PROCEEDINGS 

OF SCIENCE 



Estimate of the systematic error in tau decay due to 
duality violations 



Oscar Cata 

INFN, Laboratori Nazionali di Frascati, Via E. Fermi 40, 1-00044 Frascati, Italy 



E-mail: ocata@ifae.es 



Maarten Golterman 

Department of Physics and Astronomy, San Francisco State University, 1600 Hollow ay Ave, San 

Francis co, CA 94132, USA 

E-mail: }naarten@ stars .sfsu.edu 



Santiago Peris* 4 

Grup de Fisica Tedrica and IFAE, UAB, E-08193 Bellaterra, Barcelona, Spain 
E-mail: |peris@if ae . es 



With the help of a physically motivated ansatz, which is fitted to the data, we estimate the size of 
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estimate to be 8a s (m r ) ~ 0.003 — 0.010. Consequently, it cannot be neglected. 
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Figure 1: Left panel: Analytic structure of YU/^q 2 ) in the complex q 2 plane. The solid curve shows the 
contour used in Eq. (1.3). Right panel: Contour used in the derivation of Eq. (2.3). With the assumed 
exponential decay, Ay^q 2 ) vanish on the dashed circle with infinite radius. 



1. Introduction 

After the seminal work of Ref. hadronic tau decay has emerged as a very interesting 
physical process for studying QCD. Recent analyses find for the coupling constant a s the following 
values: 

a s (m 2 ) = 0.344 ± 0.005 exp ± 0.007 th [2] , a s (m 2 ) = 0.332 ± 0.005 exp ± 0.015th [3] , 
a s {m 2 ) = 0.321 ±0.005 exp ±0.012 th [4] , a s (m 2 ) = 0.316 ±0.003 exp ±0.005 th [5] . (1.1) 



These results show an unprecedented level of precision. It is, therefore, very important to be able 
to assess all the uncertainties involved in these analyses. 

Furthermore, Ref. after performing a very complete study of the vector and axial spectral 
functions extracted from tau decay, obtains for the gluon condensate a different number depending 
on which channel is used: 



£<GG> 



= (-0.8 ± 0.4) x 10~ 2 GeV 4 , — (GG) 

Vector K 



Axial 



(-2.2 ±0.4) x 10" 2 GeV 4 , (1.2) 



which, if taken at face value, would reflect an inconsistency in the Operator Product Expansion 
(OPE), as the gluon condensate should be the same in both channels [||]. Further determinations of 
the OPE condensates may be found in Refs. \fft. 

In the present situation we think it is time to reassess the level of precision one may expect 
from the OPE in tau decay In fact, all the analyses involve an integration of the OPE over the 
physical cut where, in fact, this expansion is expected not to converge. The central equation for 
this discussion is [[IJ, |h: 



' S0 dsP(s) -ImIV(j) = --L I dq 2 P{q 2 )n°f{q 2 ) , 



(1.3) 



where P(s) is any polynomial which may be chosen at one's convenience [ 1C], and the contour of 
integration is shown on the left panel of Fig. |T](note how the contour \q 2 \ =sq touches the physical 
cut). 
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Figure 2: Different combinations of spectral functions (V, A, and V±A) (data points in black) compared 
with Fixed-Order perturbation theory (green fiat line) and the result of our fit to Eq. (2.4) (red curve). 



Eq. ( L3) allows one to compare the OPE (on the right hand side) with the experimental data 
(on the left hand side). From this comparison, the extraction of a s and the condensates follow. 
However, if the OPE is not a good approximation on the physical cut, this extraction will have an 
error one ought to be able to estimate. Even though the polynomial P(s) may be chosen to have a 
zero on this cut, this does not guarantee that this error will vanish. This type of error is commonly 
referred to as a Duality Violation (DV). 



2. An ansatz for Duality Violations 

Assuming that the OPE is an asymptotic expansion, one may expect that the difference 



(2.1) 



will be exponentially suppressed at |<7 2 | = oo. In this case, use of Cauchy's theorem on the contour 
depicted on the right panel in Fig. [I] allows one to obtain the correction to Eq. (1.3) as 

[ S0 dsP(s) -Imll^O) = / dq 2 P(q 2 )Il° p A E (q 2 ) + 41 (*>) • (2-2) 



where [13, 14] 



poo I 

41 (*o) = - / dsP(s) -lmA VA (s) 

Jsn 71 



(2.3) 



encapsulates all the DVs. The function ImAy^ (s) in the previous equation requires knowledge of 
the spectrum up to infinite energy and is, in principle, unknown. 
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However, if indeed the OPE is an asymptotic expansion at \q 2 \ = oo , lmAy A (s) in Eq. (2.3) 
will exhibit some type of exponential falloff in s, at large values of the momentum. If, furthermore, 
we also assume that the spectrum follows some type of periodicity, such as the one found in the 
daughter trajectories in Regge theory, it becomes natural to combine an oscillatory behavior with 
the exponential falloff and assume the simple form 

1 



71 



ImAv , A (s) = k va e ^ sin (a v , A + Pv,as) 



(2.4) 



as an ansatz to fit to the experimental data. We note that a model exists [11-14] that realizes all 



these expected features of the OPE, and for which Eq. (2.4) gives the behavior at large values of s. 



3. Fitting the ansatz to tau data 



Armed with the ansatz ( J2.4[ ) we have fitted the spectral functions in tau decay in the window 1 . 1 
GeV 2 < s < m 2 .. This window has been chosen because going below the lower limit in this window 
increases the % 2 /dof very sharply, while increasing the lower limit only increases the error of the 
fit (as the data are less precise), even up to a point where the fit ceases to be meaningful. The results 
of the fits we obtained are: 



Ky = 

Yv -- 

CCv = 

Pv- 

t_ 

dof 



0.018 ±0.004 
= 0.15 ±0.15 GeV 2 
2.2±0.3 

: 2.0 ±0.1 GeV 2 
79 



K A 

7a -- 

Pa 
X 2 



dof 



0.20 ±0.06 

1.7 ±0.2 GeV" 2 

-0.4 ±0.1 

-3.0 ±0.1 GeV" 

= — ~ 0.22 . 
78 



(3.1) 



The contribution from the OPE to the spectral functions is negligible [ft]. As one can see in Fig. ||, 
and from the values of the # 2 's, the fits are very good. We remark that this is also true for the fit to 
the V + A spectral function, which shows that there is no reason why this combination of spectral 
functions should be protected from DVs, unlike what is sometimes assumed in the literature. In 
fact, because the exponent % is much larger than jy, DVs in the A channel are much suppressed 
compared to the V channel, so that V + A has the same DVs as V alone. 



Since, within the ansatz (2.4), DVs are sizeable only in the V channel, it makes sense to ask 



whether e + e~ data modify the above conclusions. Although answering this requires new assump- 
tions, the net result is that DVs tend to be decreased somewhat, but not to a level at which they 
become negligible. See the discussion in Ref. [fj]. 



4. Impact of Duality Violations on a s 



Choosing the right polynomial P(s) in Eq. (2.2) one may obtain the tau decay width on the 
left hand side of this equation as a function of the OPE parameters and the contribution from DVs 
(2.2,2.3) in terms of the fit parameters (3T). Using the total tau decay width, R r [jy, plus the 
conservative estimate for the condensates made in [ph, this leads to a shift in a s which is ^ 



8 a s (m^) ~ 0.003 - 0.010 



(4.1) 



4 



Duality Violations 



Santiago Peris 



The spread of values in (4.1) reflects the error of the fits and is insensitive to whether one uses 
the Contour-Improved or the Fixed-Order prescriptions for perturbation theory [§]. In Ref. |Q] an 
estimate of DVs was also made, with the conclusion that DVs were negligible. However, as far as 
we understand, no attempt at a detailed fit to spectral data was made there. When the fit is made, 
our conclusion is that the systematic error due to DVs is not negligible and may even get to be 
of the order of all the other systematic errors taken together. We refer to Ref. [||] for details and 
further discussions. 
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